Mapping multiple photonic qubits into and out of one solid-state atomic ensemble 
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The future challenge of quantum communication are scalable quantum networks, which require 
coherent and reversible mapping of photonic qubits onto stationary atomic systems (quantum mem- 
ories). A crucial requirement for realistic networks is the ability to efficiently store multiple qubits 
in one quantum memory. Here we demonstrate coherent and reversible mapping of 64 optical modes 
at the single photon level in the time domain onto one solid-state ensemble of rare-earth ions. Our 
light- matter interface is based on a high-bandwidth (100 MHz) atomic frequency comb, with a 
pre-determined storage time > 1 jj,s. We can then encode many qubits in short <10 ns temporal 
modes (time-bin qubits). We show the good coherence of the mapping by simultaneously storing 
and analyzing multiple time-bin qubits. 



Quantum communication [l[ offers the possibility of 
secure transmission of messages using quantum key dis- 
tribution (QKD) and teleportation of unknown quan- 
tum states [3]. Quantum communication relies on cre- 
ation, manipulation and transmission of qubits in pho- 
tonic channels. Photons have proven to be robust car- 
riers of quantum information. Yet, the transmission of 
photons through a fiber link, for instance, is inherently 
a lossy process. This leads to a probabilistic nature of 
the outcome of experiments. In large-scale quantum net- 
works [J] the possibility of synchronizing independent 
and probabilistic quantum channels will be required for 
scalability @, @. A quantum memory enables this by 
momentarily holding a photon and then releasing it when 
another part of the network is ready. In order to reach 
reasonable rates in a realistic network it will be necessary 
to use multiplexing [7j, which requires quantum memo- 
ries capable of storing many single photons in different 
modes. 

A quantum memory requires a coherent medium with 
strong coupling to a light mode. Strong and coherent in- 
teractions can been found in ensembles of atoms 0, for 
instance alkali atoms or rare-earth (RE) ions doped into 
crystals. The latter are attractive for quantum storage 
applications, as they provide solid-state systems with a 
large number of stationary atoms with excellent coher- 
ence properties. Optical coherence times of up to mil- 
liseconds 

[9] and spin coherence times > seconds [l(J have 
been demonstrated at low temperature 4 K). 

A quantum memory also requires a scheme for achiev- 
ing efficient and reversible mapping of the photonic qubit 
onto the atomic ensemble. Techniques investigated in- 
clude stopped light based on electromagnetically induced 
transparency (EIT) [Tl - 13|. Raman interactions 14- lij] 



or photon-echo based schemes (17H2l| . Much progress 



has been made in terms of quantum memory efficiency 



|13l |22| . and storage time [23j, 124J. Storage of multi- 



ple qubits is challenging, however, because it requires 
a quantum memory that can store many optical modes 
into which qubits can be encoded. Note that each pair 
of modes can encode a different qubit, or more gener- 



ally, d modes encode a qudit. 
time 0, space 



25] 



A mode can be defined in 
or frequency. Time multiplexing as 
used in classical communication has the great advantage 
of requiring only a single spatial mode (20 . [26| , hence a 
single quantum memory. Moreover, temporal modes can 
be used to define time-bin qubits [1], which are widely 
used in fiber-based quantum communication due to their 
resilience against polarization decoherence in fibers. 

This type of temporal multimode storage is difficult, 
however, due to the scaling of the number of stored modes 
Njn as a function of optical depth d of the storage medium 
20, 261. For EIT and Raman interactions N m scales as 
Vd [26J, making it very difficult to store many modes. 
Recently we proposed [20( a multimode storage scheme 
based on atomic frequency combs (AFC) with high in- 
trinsic temporal multimode capacity |20l. |26j . Using this 
method we recently demonstrated [21| that a weak co- 
herent state \oc)l with mean photon number n = |cv| 2 < 
1, can be coherently and reversibly mapped onto a YVO4 
crystal doped with neodymium ions. Later experiments 



27H29J in other RE-doped materials have improved the 
overall storage efficiency (35%) and storage time (20/xs). 
Yet, in these experiments at most 4 modes have actually 
been stored at the single photon level, thus the predicted 
[20I 26 1 high multimode capacity has yet to be shown 



experimentally. 



RESULTS 

Here we demonstrate reversible mapping of 64 tempo- 
ral modes containing weak coherent states at the single 
photon level onto one atomic ensemble in a single spatial 
mode using an AFC-based light- matter interface [20( . An 
AFC is based on a periodic modulation (with periodic- 
ity A) of the absorption profile of an inhomogeneously 
broadened optical transition \g) — > |e) (see Fig. [T|). The 
modulation should ideally consist of sharp teeth (with 
full-width at half-maximum 7) having high peak absorp- 
tion depth d, cf. Figure [TJj. Such a modulation can be 
created by optical pumping techniques (see Figure [T] and 
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FIG. 1: (a) Simplified level scheme of the Nd ions doped into Y2SiOs. We use the optical transition at 883nm between the 4 Ig 
ground state and 4 F3 excited state. The former is split into two Zeeman levels by a 0.3 Tesla magnetic field (\g) and \aux)). 
The experiment is performed on |g)-|e) where the absorption profile is shaped into an AFC by optically pumping atoms into 
\aux). The basic idea is to send in a pulse sequence on \g) - \e) that has a periodic spectral density (due its Fourier spectrum). 
Some of the excited atoms have a certain probability to spontaneously de-excite to \aux). The atoms left behind in \g) form the 
grating (see panel (b)). To build up a deep grating the sequence is repeated many times (up to the timescale of the population 
relaxation between \g) and \aux)). More details on the preparation can be found in the Supplementary Information, (b) An 
example of a generated comb with periodicity A=10 MHz. The relevant AFC parameters defined in the text are indicated, (c) 
Mapping of weak coherent states with n = 0.5 (in a single temporal mode) onto the Nd-doped crystal. Shown are two different 
experiments with A=10 MHz (dashed line offset vertically) and 1 MHz (solid line). The photons that are transmitted without 
being absorbed are detected at t=0, while the absorbed and re-emitted photons are detected around t — 1/A. The vertical 
scale has been normalized such that it yields efficiency. Inset: The overall write and read efficiency as a function of 1/A. 



Methods). This requires, however, an atomic ensemble 
with a static inhomogeneous broadening and many in- 
dependently addressable spectral channels. This can be 
found in RE-doped solids where inhomogeneous broaden- 
ing is of order 1-10 GHz and the homogeneous linewidth 
is of order 1-100 kHz when cooled < 4K. When a weak 
photonic coherent state \o)l with n < 1 is absorbed by 
the atoms in the comb, the state of the atoms can be writ- 
ten as \a) A = \G)+a\W)+0{a 2 ). Here \G) = \gi---g N ) 
represents the ground atomic state and 
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9) (1) 



represents one induced optical excitation delocalized over 
all the N atoms in the comb. In Eq. ([TJ z n is the posi- 
tion of atom n, k is the wave-number of the single-mode 
light field, 8 n the detuning of the atom with respect to 
the laser frequency and the amplitudes c„ depend on the 
frequency and on the spatial position of the particular 
atom n. The initial (at t=Q) collective strong coupling 
between the light mode and atoms is rapidly lost due 
to inhomogeneous dephasing caused by the exp(i2irS n t) 
phase factors. If we assume that the peaks are narrow 
as compared to the periodicity (i.e. a high comb finesse 
F = A/7), then S n w m n A and the W state will rephase 
after a pre-programmed time 1/A. The rephased collec- 
tive state W will cause a strong emission in the forward 
direction (as defined by the absorbed light). 



This type of photon-echo emission is also observed in 
accumulated or spectrally programmed photon echoes 
30-33], which inspired our proposal. Spectral atomic 
gratings have also been proposed [34[ and demonstrated 
[35l | for coherent optical delay of streams of strong classi- 
cal pulses. The interest in spectral gratings was recently 
renewed in the context of quantum memories, when it 
was realized how to achieve a much more efficient spec- 
tral grating than previously possible. This is possible 
due to the highly absorbing and sharp peaks in the AFC 
structure [20( . In practice the finite finesse of the comb 
still needs to be accounted for, which causes a partial 
loss of the collective state. But in Ref. [5c| we show 
theoretically that F=1Q induces a negligible loss, which 
in combination with a high optical depth d makes the 
AFC scheme very efficient. High-efficiency mapping us- 
ing high-finesse combs have been shown experimentally 
[271 . l29j . These experiments and the present work stores 
light for a pre-determined time give n by 1/A. We thus 
emphasize that we also proposed [20( and experimentally 
demonstrated [28| a way to achieve on-demand readout 
by combining AFC with spin-wave storage. On-demand 
readout is a crucial resource for applications in quantum 
networks in order to render different quantum channels 
independent. 

The multimode property of an AFC memory can easily 
be understood qualitatively. For a periodicity A and N p 
peaks, its total bandwidth is of order ~ N p A meaning 
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FIG. 2: The output from a frequency-stabilized (<100kHz) 
diode laser was split into two beams using a polarization beam 
splitter (PBS). Each beam could be amplitude, frequency and 
phase modulated using a double-pass acousto-optic modula- 
tors (AOM). One beam was used for creating the preparation 
pulses (see text), and the other one for creating the weak 
pulses to be stored (strongly attenuated using neutral density 
(ND) filters). In the weak path an additional electro-optic 
amplitude modulator (EOM) was used to create short input 
pulses for the multimode storage experiments. The paths were 
mode overlapped using a fiber-coupled beam combiner. The 
light was sent through the crystal, again in free space, in a 
double-pass setup using a Faraday rotator (FR) and a PBS. 
The output light was collected with a multimode fiber and 
detected by an Si single-photon counter (APD). Two synchro- 
nized mechanical choppers (MC) blocked the detector during 
the preparation sequence and blocked the preparation beam 
during the storage sequence, respectively. See Supplementary 
Information for more details. 



that a pulse of duration r ~ 1/ (N P A) can be stored. The 
multimode capacity stems from the fact that the grating 
can absorb a train of weak pulses before the first pulse 
is re-emitted after T = 1/A (cf. Fig. QJ). This simple 
calculation gives a multimode capacity N m oc T/t oc N p . 
Thus a comb with many peaks N p allows us to create a 
highly multimode memory in the temporal domain. In 
this context RE-doped solids are particularly interesting 
due to their high spectral channel density. 

Here we work with a neodymium-doped Y2SiC>5 crys- 
tal having a transition wavelength at 883 nm with good 
coherence properties (see Methods for the spectroscopic 
information) . This wavelength is convenient since we can 
work with a diode laser and Si based single photon coun- 
ters having low noise (300 Hz) and high efficiency (32%). 

The comb is prepared on the \g) - |e) transition by fre- 
quency selectively pumping atoms into an auxiliary state 
\aux) (see Fig. [IJ. There are different techniques for 
achieving this. For instance, by creating a large spectral 
hole, and then transferring back atoms from an auxil- 
iary state to create a comb, as used in 28j. Here we 
use a similar technique to [21(, where a series of pulses 
separated by a time r pump atoms from \g) to \aux) 



(through |e)) with a power spectrum having a periodic- 
ity 1/r = A. This technique is also frequently used in 
accumulated photon echo techniques [3(1 H^]. Here each 
pulse sequence consisted of three pulses where the cen- 
tral pulse is 7r-dephased, which has a power spectrum 
with "holes" (see Fig. [3];). A straightforward calculation 
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FIG. 3: (a) Experimental combs created using preparation 
sequences with either single (solid line) or five (dashed line) 
simultaneous pump frequencies. The frequency shifted se- 
quences allows us to enlarge the frequency range over which 
the optical pumping is efficient and thereby creating a wide 
100 MHz comb, (b) Efficiency as a function of the dura- 
tion (FWHM) of the input pulse for a single (circles) and five 
(squares) frequency preparation. As the duration decreases 
the bandwidth of the input pulse increases. The decrease in 
efficiency for short pulses is due to bandwidth mismatch for 
large bandwidths when using a single preparation frequency. 
This experiment clearly illustrates the gain in bandwidth in 
the extended preparation sequence for which only a small de- 
crease in efficiency is observed, (c) Pulse sequence for atomic 
frequency comb preparation (see text). In order to increase 
the bandwidth, the pulses are repeated with shifted frequen- 
cies. This pulse sequence was used for most of our experi- 
ments. Here it creates a comb of 100MHz bandwidth and a 
periodicity of 1MHz. The total sequence takes 16/^s and it is 
repeated around 2000 times in order to prepare the AFC. 



shows that the width of the holes in the power spectrum 
decreases with the number of pulses in the sequence. In 
this experiment three pulses were enough to reach the 
optimal comb finesse (Fss3) to achieve the maximal ef- 
ficiency for a given optical depth (see Methods). Note 
that the sum of the amplitudes of the side pulses (here 
two) should correspond to the amplitude of the central n- 
dephased pulse, in order to obtain the appropriate power 
spectrum. This rule also holds for sequences with more 
pulses. To increase the depth of the comb the sequence 
was repeated 2000 times. More details can be found in 
the Supplementary Information. 

The experimental sequence is divided into two parts: 
the preparation of the AFC (cf. above) and the stor- 
age of the weak pulses. The preparation lasts 100 ms, 
which is followed by a delay of 5 ms( ps 17Ti) to avoid 
fluorescence noise from atoms left in the excited state. 
During the storage sequence, 1000 independent trials are 
performed at a repetition rate of 200 kHz. The entire se- 
quence preparation plus storage is then repeated with a 
repetition rate of 5 Hz. An overview of the experimental 
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FIG. 4: Storage of an input state composed of 64 temporal modes. The input (left part) is a random sequence of full and empty 
time bins, where the mean photon number in the full ones is n <1. The output (right part) clearly preserves the amplitude 
information to an excellent degree. The pre-determined storage time was T=1.32 /xs, the duration of the input 1.28 {is, the 
mode separation 20 ns and the mode duration (FWHM) 5 ns. The output has been multiplied with a factor of 50 for clarity. 
The average storage and retrieval efficiency was here 1.3%. Other examples of multimode storage, e.g. with all time bins filled 
can be found in the Supplementary Information. 



set up is shown in Fig. [7j 

In Fig. QJ we show storage experiments with pre- 
determined storage times of T=100ns and 1/xs, for a sin- 
gle temporal mode. The overall in-out mapping efficien- 
cies, denned as the ratio of the output pulse counts to 
the input pulse counts, are ~ 6% and ~ 1% respectively 
(see inset of Fig. [TJ;). In the Methods section we present 
a theoretical analysis of the efficiency performance. The 
efficiency for single-mode storage is currently lower than 
has been achieved in the best-performance single-mode 
memories, e.g. 13, 16, [3, [53, 0|. But as explained 
later our interface compares very favorably to these ex- 
periments in terms of potential multimode storage effi- 
ciency. 

The main goal of the present work is to demonstrate 
high multimode storage, as theoretically predicted in 



26] . Following the discussion above, we should max- 
imize the number of peaks in the comb. This can be 
done by increasing the density of peaks in a given spec- 
tral region (i.e. increasing the storage time T) or by 
changing the width of the AFC (i.e. increasing the band- 
width). Here we fix the storage time to T—1.3 [J,s where 
we reach an efficiency of > 1% and concentrate our ef- 
forts on increasing the bandwidth. The spectral width 
of the grating is essentially given by the width of the 
power spectrum of the preparation sequence, which here 
results in a width of about 20-30 MHz. We can however 
substantially increase the total width by inserting more 
pulses in the preparation sequence, which are shifted in 
frequency (see Methods) . We thus optically pump atoms 
over a much larger frequency range. Note that the fre- 
quency shift should be a multiple of A in order to form a 
grating without discontinuities. In this way we managed 
to extend the bandwidth of the interface to 100 MHz as 
shown in Fig. [3^, without significantly affecting the AFC 
echo efficiency. This is illustrated in Figure [3Jd, where we 
show storage efficiency as a function of the duration of 
the input pulse when the preparation sequence contains 



a single or five frequencies. The maximum bandwidth 
allows us to map short <10 ns pulses into the memory. 

In addition to the present motivation for multimode 
storage, a large bandwidth is equally interesting for in- 
terfacing a memory with non-classical single-photon or 
photon pair sources. These usually have large intrinsic 
bandwidth which requires extensive filtering for match- 
ing band widths. In the present case our extended band- 
width (x5) would require a corresponding factor of less 




FIG. 5: The coherence of the multimode storage was mea- 
sured via an interference measurement, (a) The output signal 
(solid line) generated by the double- AFC scheme (see text), 
which causes an interference between consecutive modes. The 
input sequence (not shown) is a series of weak coherent states 
(n sal.8) of equal amplitude (c| = c^) where the relative 
phase difference between modes <j>k+i — <j>k ranges from -7r/6 
to 8n/6 with a step of n/6. This allows us to capture a com- 
plete interference fringe in one measurement. It also clearly 
demonstrates the preservation of coherence over the complete 
multimode output. By changing the detuning of the centre of 
the second AFC with respect to the carrier frequency of the 
light, we can impose an additional relative tt phase [20I ] on 
the corresponding output (see Supplementary Information). 
This shifts the interference fringe with half a period (dashed 
line), (b) The corresponding net interference visibilities are 
86%±3 (open circles) and 85%±3% (filled squares), with de- 
tector noise (dashed line) subtracted. The uncorrected raw 
visibilities were 78%±3% and 76%±2%. 
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filtering. 

We show the high multimode capacity of our interface 
by storing 64 temporal modes during a pre-determined 
time of 1.3 microseconds (see Fig. [4]), with an overall 
efficiency of 1.3%. This capacity is more than an order 
of magnitude higher than previously achieved for mul- 
tiplexing a quantum memory in a single spatial mode 
16l 21 1 . As shown we can store a random sequence of 
weak coherent states. Storage of random trains of single 
photon states has been proposed for multiplexing long- 
distance quantum communication systems based on so 
called quantum repeaters 0, @] . The maximum rate of 
communication would then be proportional to the num- 
ber of modes that can be stored Q. Our experiment 
clearly shows the gain that can be made using an AFC- 
based quantum memory. It thus opens up a route to- 
wards achieving efficient quantum communication using 
quantum repeaters. 

It is now possible to use consecutive temporal modes, 
e.g. modes \k) and \k + 1), to encode time- bin qubits 



c k \k) + c 



,i<t>kn 



| to), in which case a good coherence be- 



tween modes is crucial. The coherence can be measured 
by preparing superposition states and performing pro- 
jective measurements using an interferometric set up. 
Projective measurements on time-bin qubits is usually 
performed using an unbalanced Mach-Zehnder interfer- 
ometer (MZI) where consecutive time-bin are interfering 
[l[ . We can perform the same task with our light- matter 
interface by using a double-AFC scheme (with Ai and 
A2) as shown in [2lj. In short, the difference in de- 
lay 1/Ai — I/A2 plays the role of the delay in a un- 
balanced MZI. We observe excellent coherence over all 
modes with an average visibility of V = 86% ± 3%, see 
Fig. corresponding to a conditional qubit fidelity of 
F = (l + V)/2«93%. 

To further illustrate our ability to store multimode 
light states we create a light pulse with a random am- 
plitude modulation. As shown in Fig. [6]we can faithfully 
store this kind of light pulses. The possibility of stor- 
ing weak arbitrary light states using photon-echo based 
schemes was pointed out already by Kraus et al. [18| ■ We 
believe that this work, where complex phase and ampli- 
tude information are reversible and coherently mapped 
onto one atomic ensemble, is the first experimental re- 
alization showing these properties at the single photon 
level. 



DISCUSSION 

For multimode storage the efficiency of our interface 
would outperform the current EIT and Raman based 
quantum memories in homogeneously broadened media, 
although impressive efficiencies have been achieved for 
single-mode storage [8, 13, 22 1. This is due to the poor 
scaling of the efficiency as a function of the number of 
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FIG. 6: Mapping of a 1 /is long input pulse with randomly 
varying amplitude. As seen the overlap between the normal- 
ized input (dashed line) and output (solid line) pulses is ex- 
cellent. The total average number of photons in the input 
pulse is n « 4 



modes for a given optical depth 26] . It also compares fa- 



vorably to the recent few modes storage experiment [16 1 
using the gradient echo memory (GEM), another echo 
based storage scheme, also due to the scalin g of mode 
capacity for a given optical depth (N m ~ d) fj\ |26| . Still, 
an increase in storage efficiency and on-demand read out 
is necessary for applications in quantum communication. 

The next grand challenge is to combine multimode 
storage, high efficiency [29[ and on-demand read out [28[ 
in one experiment. The immediate efforts will most prob- 
ably by devoted to praseodymium and europium-doped 
Y2SiC>5 crystals, where the ground state manifold has 
the necessary number of spin levels (three levels) for im- 
plementing the on-demand readout. The recent achieve- 
ments in Pr-doped Y2SiOs crystals are very encouraging 
[28l [29|, although the bandwidth was limited to a few 
MHz due to the hyperfin level splitting. Europium-doped 
Y2SiC>5 has the potential of offering higher band widths 
(up to 70 MHz) and narrower comb peaks, which results 
in higher multimode capacity [2(J. In order to exploit 
the high-bandwidth results reported in this work, using 
neodymium-doped crystals, one needs to find a third spin 
level with a long spin coherence lifetime. An interesting 
path forward is to investigate neodymium isotopes with 
a hyperfine structure ( 143 Nd and 145 Nd) (36|. Recent 
results on a similar system 37], 167 Er 3+ :CaW04, show 
coherence times approaching 100 fxs for hyperfine transi- 
tions. Clearly this path requires extensive spectroscopic 
studies in order to optimize the spin population and co- 
herence lifetimes. But it is very interesting since it opens 
up several material candidates (e.g. doped with Erbium 
[381 ] and Neodymium) for quantum memory applications. 

To summarize we have demonstrated the reversible 
mapping of up to 64 optical temporal modes at the single 
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photon level onto one solid state atomic ensemble. We 
have demonstrated that the quantum coherence of the 
stored modes is preserved to a high extent. The differ- 
ent modes can then be used to encode multiple time-bin 
photonic qubits. Alternatively, they could also be con- 
sidered as high-dimensional qudits states. This opens up 
possibilities to store higher dimensional quantum states 
such as entangled qudits encoded in time bin bases. Our 
experiment opens the way to multi-qubit quantum mem- 
ories, which are a crucial requirement for realistic quan- 
tum networks. 



METHODS 

Sample 

The sample is a 10mm long neodymium-doped yttrium orthosil- 
icate crystal(Nd 3_, ~ :Y2SiOs) with a low Nd 3+ concentration of 30 
ppm. The inhomogeneous broadening of the 4 Is>- 4 F3 absorption 

2 2 

line is around 6GHz and the optical depth 1.5 for this sample. By 
using a double-pass set up through the crystal we could increase 
the optical depth to 3. We measured an excited state lifetime of 
Ti=300^ts using fluorescence spectroscopy and stimulated photon 
echoes. With conventional photon echoes (two-pulse) we measure 
a homogeneous linewidth of 3.4kHz (T2 = 92.7/^s). Each level is 
a Kramer's doublet which split into two spin states in a magnetic 
field. For the field orientation used in this experiment we measured 
g factors of g g = 2.6 and g e = 0.5. In a 300mT magnetic field, the 
excited states were separated with 2GHz. We measured a ground 
state Zecman population relaxation lifetime, by spectral hole burn- 
ing (SHB), of around Ti^ = 100ms. In the SHB measurements we 
also observed a superhyperfine interaction of Nd ions with yttrium. 
This causes additional spectral side holes at around 640kHz (for the 
present magnetic field), thus the effective homogeneous linewidth 
is around 1 MHz. This was our main limitation for the efficiency 
of our light-matter interface since it affected our ability to create a 
good comb for the longer storage times (1/A lfis). 
Storage efficiency analysis 

The efficiency can be calculated theoretically using the formula 
[UHl V ~ W 'F) 2 e~ d l ' F 'e~~> 'I ' p2 ' e~ d o . The different terms can be 
given a qualitative understanding. The first term represents the 
collective coupling, the second the re-absorption of the re-emitted 
light, the third is an intrinsic dephasing factor due to the finesse 
and the last term a loss due to an absorption background do- For 
the comb with A=10MHz we measure d £3 1.7,-F 2.7 and do & 0.5 
(see Fig. [lb), resulting in a theoretical efficiency of r\ m 5% in close 
agreement with the experiment (see Fig. [T};). The major limiting 
factor here is do (caused by imperfect preparation of the comb) 
and then the optical depth of the comb d (the finesse being close 
to optimum for this d |20|0 . The decrease in efficiency for longer 
storage times (see inset of Fig. [T};) is principally due to an increase 
in the background absorption do and an accompanying decrease 
in the peak absorption d. This in turn is caused by the effective 
spectral resolution of 1 MHz in the optical pumping, which is a 
limitation of the present material (see Methods above). However, 
the storage efficiency is between one and two orders of magnitude 
higher than what we achieved in the material Nd:YVC>4 |2ll . which 
we attribute to an improvement in optical pumping in this Nd- 
doped material. 
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Supplementary Information 

We provide additional information for our manuscript. This 
material includes a detailed description of the experimental setup 
and of the atomic frequency comb preparation, as well as additional 
results. 



Experimental details 
Setup 

The experimental setup is shown in Fig[7] The light is emitted 
by an external cavity diode laser (Toptica) at 883nm and is then 
split with a polarizing beam splitter (PBS) into two paths. In 
both paths, the amplitude, frequency and phase of the light can 
be modulated using acousto-optic modulators (AOM, AA-Opto- 
Electronics) in double pass configuration. The path on the top is 
used to prepare an AFC into the sample. To achieve this, the AOM 
creates a series of pulses with a characteristic spectrum explained in 
the next section. The other path is used to create the weak pulses 
of light at the single photon level that are going to be mapped 
onto the crystal. With an AOM and neutral density(ND) filters 
we create pulses containing less than one photon in average. The 
AOMs have a good extinction ratio (> 40 dB) and allows us to shift 
the frequency or the phase of the light. However, to create very 
short pulses(few ns), we use an additional a 20 GHz electro-optical 
modulator (EOM, EO Space). 




FIG. 7: Experimental setup : see text for details. The fol- 
lowing abbreviations are used : polarizing beam splitter (PBS), 
acousto-optic modulator (AOM), neutral density filter(ND filter), 
mechanical chopper (MC), polarization controller (PC), electro- 
optical modulator (EOM), Faraday rotator (FR) and avalanche 
photo-diode (APD). 

Both path are then combined using a 90/10 fiber coupler and 
sent into free space into the sample cooled at 3K using a pulse tube 
refrigerator (Oxford Instruments). The PBS combined to a A/2 
wave-plate allow us to adjust the polarization of the light parallel to 
the Dl crystallographic axis (optical axis). In order to increase the 
available optical depth, we implemented a double-pass setup using 
a Faraday Rotator (FR). In this way, we have a configuration where 
the polarization of the light remains constant when propagating in 
double pass through the crystal, but is rotated by 90 degrees after 
the second pass in the FR, such that the the light emitted by the 
crystal can be separated from the input light, thanks to a PBS. 

The light is then coupled to a multimode optical fiber and de- 
tected by a single photon Si Avalanche Photo-Diode (APD) with 
32% efficiency. The mechanical chopper(MC) in the preparation 
path is used to block the leakage of the preparation AOM when 
the APD is detecting echoes. The second MC is used to protect 
the APD when strong preparation pulses are used. The transmis- 
sion between the input of the cryostat and the APD was typically 
between 25% and 30 %. Finally, a third path (partially shown), is 
used to actively frequency stabilize the laser to less than 100kHz 
using the Pound-Drever-Hall technique. The frequency reference is 



given by a spectral hole in the crystal. The stabilization light is 
sent into the same sample, but with a slightly different angle. 

The experiment is repeated every 200ms, allowing sufficient time 
for the ground state population to relax between experiments. The 
first 100ms were used for the preparation sequence to create an 
AFC. Then, we waited 5ms to avoid fluorescence during the stor- 
age and retrieval sequence, due to atoms left in the excited state 
after the preparation of the AFC. After that, we perform ~ 1000 
independent storage trials separated by 5fis. Each trial contains 
the particular pulse sequence to be stored and the re-emitted AFC 
echoes. Depending on the mean photon number per pulses, the to- 
tal integration time to accumulate sufficient statistics varied from 
30s to lOmin. 



Comb Preparation 

We now explain in more detail the preparation sequence allowing 
us to create a desired AFC. Similarly to [y, we send series of pulses 
separated by a time r, which create a periodicity of A = 1/r in the 
absorption profile. By sending repeatedly such sequence, this mod- 
ulation will take the form of a comb with sharp peaks. However, 
instead of sending only pairs of pulses as in [lj, we extended the 
method by allowing the possibility of sending N pulses. The spec- 
trum of this series of N pulses is a a comb of periodicity A = 1/r. 
However, we want to remove from the initial state only atoms that 
are not in the desired comb. The spectrum of the sequence must 
thus be the inverse of a comb : a series of holes separated by A. 
To achieve this, the central pulse must be 7r-dephased, and with a 
pulse area equal to the sum of the other pulses (see Fig. 2). 

The property of the comb can be determined by the characteris- 
tic of the pulses sequence. As already mentioned, the periodicity A 
is given by the time r between pulses. The duration of one pulse(or 
its spectrum) will determine the bandwidth of the whole comb. Fi- 
nally, the width of each peak will be inversely proportional to the 
duration of the whole pulse sequence. Similarly, we can say that 
the finesse of the comb will increase with the number of pulses (F~ 
N pu i ses ). Note that these rules are true for the spectrum of the 
light, but the process of spectral hole burning being more complex, 
the property of the AFC can be different. For example, the width of 
each peak is limited by the crystal properties, such as homogeneous 
linewidth and superhyperfine interaction. 

The comb preparation sequence for most of our experiments is 
shown in Fig|8] It takes 16^ts and it is repeated around 2000 times. 
The available optical depth being relatively small (d=3), the op- 
timal finesse to get the maximum storage end retrieval efficiency 
was about 3 Q . It was thus sufficient to send only series of three 
pulses (N=3)to create the desired comb. The spectrum of the comb 
created in this way was around 20MHz. It was possible to ex- 
tend this bandwidth by sending other series of pulses at a shifted 
frequency(±20MHz,±40MHz,etc. . . ). Since there is no interference 
between different frequencies, the pulses can be sent within the co- 
herence time Ta. This means that the whole preparation sequence 
take the same time as with a single frequency, and does not require 
more Rabi frequency. 



Double read-out using double AFC 

For the interference experiment, the stored excitations must be 
read-out at different times. For that purpose, we need to create 
multiple AFCs, with different periodicity A; and possibly different 
phases. This can be realized by using a preparation pulse sequences 
with different pulse separations. Suppose that we want to create 
two AFCs with periodicity Ai and A2. The necessary sequence is 
shown in Fig|9] We send a pulse sequence with N pulses separated 
by n (with the central pulse at t=0), superposed with N pulses 
separated by T2 (with the central pulse also at t=0). Thus, the total 
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FIG. 8: Pulse sequence for atomic frequency comb preparation. 
In order to increase the bandwidth, the pulses are repeated with 
shifted frequencies. This pulse sequence was used for most of our 
experiments. Here it creates a comb of 100MHz bandwidth and 
a periodicity of 1MHz. The total sequence takes 16^ts and it is 
repeated around 2000 times in order to prepare the AFC. 



with 1MHz peak separation (corresponding to 1 fis storage time) 

and 20 MHz and 100 MHz bandwidth, respectively. As noted, be- 
fore the enlargement of the bandwidth does not affect the center of 
the comb. Thus it allows us to increase the number of mode for a 
fixed storage time without affecting the efficiency (see Fig. 2 of the 
manuscript). However, compared to combs with bigger peak sepa- 
ration(shorter storage time), the height of the peaks has decreased, 
and the absorption background has increased, which explained the 
decay of efficiency with storage time. As explained in the paper, 
the number of mode is proportional to the number of peaks in the 
comb. Here we create more than 100 peaks which illustrate the 
multimode performance of our experiments. 



number of pulses will be 2N-1, and the pulse area of the central one 
is equal to the sum of the others. When we send a storage pulse 
in the sample, two echo will be emitted, after a time t\ and tq.. 
In our case, we use N=3, such that the total number of pulses for 
creating two AFCs with different periodicity is 5 (for one pumping 
frequency). Note that when we superpose two combs, some peaks 
of the two AFCs can be at the same position and be summed. This 
means that some peaks will be higher than others. Thus, using 
this method the two AFCs cannot be created independently (i.e. 
with two independent pulses series) or the absorption profile will 
not correspond to the sum of two AFCs and we will face additional 
echoes. So we must create directly two AFCs in one sequence. 

Finally, we would like to be able to induce a phase rf> in one of the 
AFCs, in other world to shift the frequency of the comb Q ■ To do 
so, we must add a phase in each pulse preparing the corresponding 
comb (Sec Fig. 3). If we label the pulses on the right of the central 
one with k=l,2,3. . . and the pulses on the left with k=-l,-2,-3. . . , 
then the phase in the pulse k must be k</>. 
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FIG. 10: Experimental atomic frequency combs with 1 MHz peak 
separation for single(green line) or five(black line) different simul- 
taneous pump frequencies. 



Finally, we show two more example of storage of 64 temporal 
modes with different inputs. In Fig llll all the modes are full, 
which allows us in principle to store 32 time-bin qubits. In Fig |12l 
we modulated the amplitude of the input pulses. We note that 
even if the input mode changes, we do not need to adapt the AFC 
preparation. Indeed, the efficiency is constant for all modes, and 
the echo always closely follow the input pulses. 



FIG. 9: Pulse sequence for the preparation of two atomic frequency 
combs, as required for the double read out for the interference 
experiments. In this example, the two AFCs have a periodicity 
1/ti and I/V2. The first AFC will also get a phase cj>. 



Additional results 



We present here some additional results that illustrate the mul- 
timode capacity of AFC. In Fig. 1101 we show experimental combs 
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FIG. 11: Storage of 64 consecutive pulses. An efficiency of 1.4% was measured. The mean photon number per pulse is n ^ 0.5 
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FIG. 12: Storage of 64 pulses with a modulating amplitude. An efficiency of 1.6% was measured. The mean photon number in the biggest 
pulses is n fs 0.9 



